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ABSTRACT 

Recent high-quality photometry of many star clusters in the Magellanic Clouds with ages of 
1-2 Gyr revealed main sequence turnoffs (MSTOs) that are significantly wider than can be 
accounted for by a simple stellar population (SSP). Such extended MSTOs (eMSTOs) are 
often interpreted in terms of an age spread of several 10^ yr, challenging the traditional view 
of star clusters as being formed in a single star formation episode. Li et al. and Bastian & 
Niederhofer recently investigated the sub-giant branches (SGBs) of NGC 1651, NGC 1806, 
and NGC 1846, three star clusters in the Large Magellanic Cloud (LMC) that exhibit an 
eMSTO. They argued that the SGB of these star clusters can be explained only by a SSR 
We study these and two other similar star clusters in the LMC, using extensive simulations 
of SSPs including unresolved binaries. We find that the shapes of the cross-SGB profiles of 
all star clusters in our sample are in fact consistent with their cross-MSTO profiles when 
the latter are interpreted as age distributions. Conversely, SGB morphologies of star clusters 
with eMSTOs are found to be inconsistent with those of simulated SSPs. Finally, we create 
PARSEC isochrones from tracks featuring a grid of convective overshoot levels and a very 
fine grid of stellar masses. A comparison of the observed photometry with these isochrones 
shows that the morphology of the red clump (RC) of such star clusters is also consistent with 
that implied by their MSTO in the age spread scenario. We conclude that the SGB and RC 
morphologies of star clusters featuring eMSTOs are consistent with the scenario in which the 
eMSTOs are caused by a distribution of stellar ages. 

Key words: galaxies: Magellanic Clouds, star clusters; globular clusters: general; stars: 
Hertzsprung-Russell and colour-magnitude diagrams 


1 INTRODUCTION 

Until recently, globular clusters (GCs) were thought to be sim¬ 
ple stellar populations (hereafter SSPs), consisting of thousands 
to millions of coeval stars with the same chemical composition. 
In the past decade, however, a consensus has emerged that GCs 
typically harbor multiple stel lar populations featuring several un¬ 
expected characteri stics (e.g.. [Bedin et alJ[2004l : IPiotto et alj|2007l : 
iGratton et al.l l20lT and references therein). Multiple sequences 
in several major features of colour-magnitude diagrams (CMDs) 
are now commonplace in massive Galactic star clusters, e.g, in 
their main sequence (MS), sub-giant branch (SGB), and red giant 
branch (RGB). These multiple sequences are especially evident in 
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Hubble Space Telescope (HST) photometry using colours involv¬ 
ing passbands that cover molecular features of OH, NH, and CN 
in t he near-ultraviolet part of the electromagnetic spectrum (see , 
e.g.JPiotto et alj2012l:lBellini et alj20lllMilone et alj20ll|2013l. 

|2()li|D otter et alj2015l). 

Meanwhile, recent spectroscopic surveys established that light 
elements such as N, O, and Na show large star-to-star abundance 
variations (often dubbed “Na-0 anticorrelations”) within virtually 
all G alactic GCs studied to date in sufficient detail jCarretta et al.l 
I 2 OIOI and references therein). The chemical processes involved in 
causing the light-element abundance variations have largely been 
identified as proton capture reactions at high temperature (T ^ 
2 X 10 ’ K), such as the CNO and NeNa cycles. Currently, the lead¬ 
ing candidates for polluting sources (“polluters”) are stars in which 
such reactions occur readily and which feature slow stellar winds 
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so that their ejecta are relatively easy to retain within the poten¬ 
tial well of massive clusters: (1) intermediate-mass asymptotic gi- 
ant branch (AGB) and sup er-AGB stars (4 < M-jM^ < 10; e.g., 
IP’Antona & Venturall2007h, (2) rapidly rotating m assive stars (of¬ 
ten referred to as “FRMS”; Decressin et alJl2Q07h and (3) massive 
binary stars l lde Mink et al.ll2009t) . 

In the two currently favored formation scenarios, these chem¬ 
ical anticorrelations are due to stars having either formed from 
or polluted by gas that is a mixture of pristine material and ma¬ 
terial shed by suc h polluters. In the “extended star formation” 
scenario (see, e.g., Ip Ercole et alj|2008l ; IConrov & Snergellboill ; 
IVancarce & Catelanll20Ilh ~ the abundance variations are caused by 
a second generation of stars that formed out of gas clouds that were 
polluted by winds of first-generation stars to varying extents, dur¬ 
ing a period spanning up to a few hundreds of Myr, depending 
on the nature of the polluters and the depth o f the potential well . 
The alternative “early disc accretion” scenario jBastian et alj2013h 
does not involve extended star formation: the polluted gas is instead 
produced by FRMS and massive binary stars, and accreted by low- 
mass pre-main-sequence stars during the first « 20 Myr after the 
formation of the star cluster. 

In the context of the nature of Na-O anticorrelations in 
Galactic GCs, the recent discovery of extended main sequence 
turn-offs (hereafter eMSTOs) in interm ediate-age (1-2 Gyr old) 
star c l usters in the Magell a nic Clouds (jMacte^^_Brob^Nielsen| 
20071: iMackev et al.l |2^: iGlatt et al.l l2008l: iMilone et al.l l2009l : 


Goudfrooij et al.ll200^ ' 2011al lbh has generated much interest. Sev¬ 


eral investigations suggested that the eMSTOs are due to the pres¬ 
ence of multiple stellar populations span ning an age interval of sev¬ 
eral 10^ yr within these clusters (see also Rubele,.Kgrber,_^(jirardi| 


201(]l:lRubele et al.ll201ll:lKeller. Mackey. & Da Costall201 ll |2012|: 


Girardi et alj 2013£J : Correnti et al.l 20141 : Goudfrooii et al.ll2014 ). 

The leading alternative theory for the cause of the eMSTO 
phenomenon is that it is due to spreads in rotation veloc¬ 
ity a mong turnoff stars (hereafter the “stel l ar rotation” sce- 
nario: Bastian_^deMink l2009l : llJet_alJ^ 20121; Yang et al]l2013l ; 
iLi. de Griis. & Dendl2014L but see lGirardi et “^12011 . 

Due to the important implications of the answer to the ques¬ 
tion whether or not eMSTOs in intermediate-age clusters are due 
to a significant range in stellar ages, the issue is being studied from 
a variety of angles. Recent investigations looked for evidence of 
extended star formation (or a significant lack thereof) in young 
massive star clusters in the Large Magellanic Cloud (LMC), with 
see mingly conflicting re sults Isee lSastian & Silva-Villa|[2013l ver¬ 
sus jConenti^^r^lTOl^)- Others studied relations between prop¬ 
erties of eMSTOs and dynamical properties of the clusters, find¬ 
ing a correlation between eMSTO width and cluster escape ve¬ 
locity which seems to be understood most str aightforwardly in 
terms of the extended star formation scenario (iGoudfrooii et al.l 
[ 20 TI hereafter G-h 14). Howev er, result s with the opposite conclu- 
sion were recently repo rted bv iLi et ^ (2014. hereafter L-fl4) and 
ISastian & NiederhofeJ (2015, hereafter BN15). L-l-14 studied the 
SGB morphology of NGC 1651, an eMSTO cluster of age ~ 2 Gyr, 
and concluded that it can be explained only by a single-age stellar 
population. Showing results of new calculations of stellar tracks 
in the MSTO-SGB region of the CMD with vs. without stellar ro¬ 
tation, L-l-14 argued that eMSTOs in intermediate-age clusters are 
most likely caused by a range of stellar rotation velocities within 
such clusters. BN15 performed a similar study of the SGB mor¬ 
phologies in NGC 1806 and NGC 1846, two other eMSTO clus¬ 
ters, again concluding that age spreads are unlikely to be the cause 
of the eMSTO phenomenon. 


In the current paper, we investigate the claims of L-l-14 and 
BN15 by means of an independent study of the SGB properties of 
five intermediate-age star clusters in the LMC (including the clus¬ 
ters studied by L-l-14 and BN15). After a brief description of the 
data and SSP models used in this paper in Sect.|2 we discuss prop¬ 
erties of cross-SGB magnitude distributions for these star clusters 
and highlight differences with the methodologies used by L-l-14 and 
BN15 in Sect. 13 Sect, [previews the relevance of convective over¬ 
shoot to the morphology of the CMD of intermediate-age clusters, 
and discusses our cross-SGB distributions when compared with a 
set of isochrones in which the implementation of the dependence of 
convective overshoot on stellar mass is different from that used in 
the isochrones used by L-l-14 and BN15. Sect.|4]also includes an in¬ 
vestigation of the claim by BN 15 that the morphologies of the red 
clumps in these star clusters are inconsistent with the age spread 
scenario. Our conclusions and their implications are summarized 
in Sect. 13 


2 DATA AND MODELS 


We use the HST photometry of inte rmediate-age star clusters in the 
LMC that was described i n detail in Goudfrooii et aljbOllaL here¬ 
after G-l-lla), G-l-14, and ICorrenti et alj | |2014|) . Briefly, the data 
from G-l-lla used here involves multi-exposure photometry taken 
with the Wide Field Camera (WFC) of the Advanced Camera for 
Surveys (ACS) using the F435W and F814W filters {HST program 
10595, PI: P. Goudfrooij). The data from G-l-14 used here involves 
similar photometry, now taken with the Wide Field Camera #3 
(WFC3) using the F475W and F81 4W filters (HST progr am 12257, 
PI: L. Girardi), while the data from lCorrenti et al.l 1 I2OI4I) used here 
involves ACSAVFC photometry using the F555W and F814W fil¬ 
ters with one exposure per filter {HST “snapshot” program 9891, 
PI: G. Gilmore). 

From this sample we first select the star clusters in the LMC 
for which the eMSTO feature is clearly established. Specifically, 
we select star clusters for which G-l-14 measured full-width at half 
maximum (FWHM) values of their MSTO cross-cut of at least 350 
Myr when expressed as an age range (i.e., their “pseudo-age dis¬ 
tributions”). Since we are interested in measuring the SGB mor¬ 
phologies of such clusters, we subselect clusters for which their 
CMD includes at least 35 stars in the SGB regiorfl This selection 
procedure yields the three massive clusters NGC 1783, NGC 1806, 
and NGC 1846. To this sample we add NGC 1651 (which also fea¬ 
tures an eMSTO), the cluster for which lETet alj | |201^ claimed that 
its SGB morphology can be explained only by a single-age stellar 
population. Finally, we add IC 2146, a cluster whose mass and age 
are very similar to those of NGC 1651 although it does not show 
evidence for an eMSTO. We find that the SGB properties of IC 
2146 reveal information that is relevant to the question whether or 
not (and, if so, how) SGB morphologies of eMSTO clusters can be 
used to rule out extended star formation periods. Relevant proper¬ 
ties of the star clusters in our sample are listed in Table[T] 

To avoid issues related to the presence of a significant number 
of LMC field stars in relevant parts of the CMDs, we only consider 
stars within the core radius of the clusters. The only exception to 
this is the case of NGC 1651 for which the combination of cluster 
mass and core radius would yield too few stars in the SGB region 


' We found that clusters with fewer than 30-35 SGB stars showed signifi¬ 
cant stochastic fluctuations in their cross-SGB magnitude distributions (see 

Sect.|3). 
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Figure 1. Illustration of the assignment of relative masses to individual SSP ages (i.e., isochrones) for the “multi-SSP” simulations of the 4 star clusters in our 
sample that feature eMSTOs (see Sect. Panel (a): The black solid line represents the pseudo-age distribution of NGC 1651 taken from G+14, expressed 
as probability density function (PDF). The histogram in light grey shows the relative weights of the different age bins used during the simulations. Panel (b): 
Same as Panel (a), but now for NGC 1783, using its pseudo-age distribution taken from G+lla, using ages from the M-i-08 isochr ones. Panel (c): Same a s 
Panel (b), but now for NGC 1806. Panel (d): Same as Panel (b), but now for NGC 1846, using its pseudo-age distribution taken from lGoudfrooii etalJtoOgl) . 


Table 1. Properties of stai* clusters in our sample. 


Name 

(1) 

log Aid 
(2) 

rc 

(3) 

^eff 

(4) 

Ref. 

(5) 

Age Ref. 
(6) (7) 

NGC 1651 

4.91 ±0.06 

4.57 ±0.36 

12.82 ±2.01 

1 

2.00 

1 

NGC 1783 

5.42 ±0.11 

10.50 ±0.49 

11.40 ±2.24 

1 

1.70 

2 

NGC 1806 

5.10±0.06 

5.91 ±0.27 

9.04 ±1.24 

1 

1.65 

2 

NGC 1846 

5.24 ±0.09 

8.02 ±0.49 

8.82 ±0.68 

1 

1.65 

3 

IC2146 

4.49 ±0.07 

8.89± 1.36 

12.53 ±1.92 

4 

1.90 

4 


Notes. Columns: (1) Name of star cluster. (2) Logarithm of cluster mass 
in Mq. (3) Core radius in pc. (4) Effective radius in pc. (5) Reference 
for cl uste r mass and radii (1 = G+14, 2 = G+lla, 3 = iGoudfrooii et all 
I 2 OO 9 I 4 = ICorrenti et ^l2014h . (6) Age in Gyr used for the “single-SSP” 
simulations in Sect.[2 (7) Reference for age data. 


to produce statistically robust results. Hence we use all stars within 
the effective radius for NGC 1651. The influence of field stars in the 
SGB region is negligible in all cases, as can be appreciated from the 
various CMD plots (panels (a) in Fi gs.[^and[4l-[^. 

We adopt the isochrones of iMarigo et alj j2008l . hereafter 
M+08) for much of our analysis in this pape r, since these 
isochr ones were also used by G+lla, G+14, L+14. I^rrenti et ^ 
( l2014h . and BN15. For each cluster, we adopt the values for Z (al¬ 
ways 0.008), distance (m-M)o, best-fittin g age, and foregroun d 
reddening (Av) found by G+lla, G+14, and lCorrenti et al.l d201 4h. 
The population properties of these clusters were determined from 
the brightnesses and colors of the MSTO and the RGB bump, and 
the slope of the RGB. The SGB morphology was not involved in 
the determination of cluster ages. 

3 SGB MORPHOLOGIES 
3.1 Method 

To test whether the SGB morphologies of the clusters in this paper 
are best described by a single-age SSP or by a distribution of ages 
similar to those indicated by their MSTO morphologies (i.e., the 
pseudo-age distributions shown in G+lla and G+14), we create 
cross-SGB magnitude distributions (hereafter called “cross-SGB 


profiles”) of the observations and compare them with such profiles 
derived from a set of Monte Carlo simulations of SSPs (including 
unresolved binary stars) as described below. 

Candidate SGB stars are selected by means of a parallelogram 
in the CMD, illustrated in panels (a)-(c) of Fig.j^by black dashed 
lines. The blue edge of the parallelogram is placed at a colour that 
is clearly redder than the MSTO feature, while the red edge of the 
parallelogram is placed at the colour reached by the SGB stars with 
the minimum luminosity in the relevant isochrones. The placement 
of the top and bottom of the parallelogram is guided by the distri¬ 
bution of stars in the various SSP simulations, making sure all ar¬ 
tificial SGB stars are captured (including unresolved binaries) and 
verifying that the cluster’s SGB stars are captured as well. Cross- 
SGB profiles are then created by measuring every star’s magnitude 
offset from the isochrone wi th the best-fitting (me an) age as estab¬ 
lished in G+lla, G+14, and ICorrenti et al.l l l2014h at the colour of 
the star in question, using linear interpolation in colour space. 

SSP simulations are co nducted by popu lating isochrones with 
stars randomly drawn from a lSalpetea d 1 9550 IMF between the min¬ 
imum and maximum stellar masses in the isochrone. We add unre¬ 
solved binary companions to a fracti on of the stars using th e binary 
fractions found by G+lla, G+14, and ICorrenti et al.l ( l2014h . in con¬ 
junction with a flat primary-to-secondary mass ratio distribution. 
Finally, we add random photometric errors whose dependence on 
the location in the CMD was d erived from th e artifi cial star tests 
described by G+lla, G+14, and ICorrenti et alj ( |2014|) . 

We perform two distinct sets of SSP simulations for each clus¬ 
ter: (1) a set of “single SSP” simulations, using the isochrone with 
the best-fitting (mean) age; (2) a set of “multi-SSP” (or “composite 
SSP”) simulations. These use the isochrones whose ages encom¬ 
pass the range of ages indicated by the clusters’ pseudo-age distri¬ 
butions that were derived from MSTO crosscut profiles by G+lla 
and G+14. We employ an age resolution of 0.1 Gyr. For each age, 
the relative number of stars used in the simulation is set by the me¬ 
dian amplitude of the pseudo-age distribution for that age bin, in 
probability density units (this procedure is illustrated in Figure fTpl 
The level to which such “multi-SSP” simulations approximate the 
relevant CMD features of the star cluster in question is illustrated 


^ This procedure is somewhat different from that used by BN15, who fitted 
single gaussians to the clusters’ binned MSTO crosscut profiles. 
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Figure 2. Panel (a): the CMD of NGC 1651 in the MSTO-SGB region (taken from G+14) along with lMarigo et alj i2QQ8l) isochrones for Z = 0.008 and ages 
1.7, 1.8, 1.9, 2.0, and 2.1 Gyr. Black dots represent stars within the effective radius of NGC 1651. Red dots indicate stars in the outskirts of the HST image 
with a total area equal to that within the effective radius of NGC 1651. The SGB selection box is outlined by black dashed lines. For reference, the MSTO 
crosscut parallelogram used by G+14 to measure the pseudo-age distribution of NGC 1651 is shown in magenta. Panel (by. CMD of a single-SSP simulation 
of NGC 1651 using the 2.0 Gyr isochrone (see Sect. 1 ^ . Single stars are shown as blue dots while unresolved binaries are shown as orange dots. The black 
solid line within the SGB selection box in panel (b) represents the SGB portion of the isochrone for the best-fitting (mean) age. The black arrows in panel (b) 
indicate the direction of the positive X axis in the cross-SGB profiles shown in Figures[3-[2l Panel (c}\ similar to panel (b) except that the dots now represent 
a “multi-SSP” simulation of NGC 1651, described in Sect. [XT] Simulated stars with ages of 1.7, 1.8, 1.9, 2.0, 2.1, and 2.2 Gyr are shown in green, blue, red, 
black, purple, and cyan, respectively. Panel (d): the black solid line represents the pseudo-age disti'ibution of NGC 1651 as measured from the MSTO crosscut 
by G+14 (cf. Fig la), expressed as a normalised density function, while the yellow region indicates its 68% confidence interval. For comparison, the red dashed 
line represents the corresponding pseudo-age distribution of the “multi-SSP” simulation described in Sect. IXTIforNGC 1651. 


in Fig.l^for the case of NGC 1651. In particular, panel (d) of Fig. 
shows a comparison of the pseudo-age distribution of NGC 1651 
as derived from its MSTO crosscut profile by G-l-14 (i.e., the curve 
shown also in Fig. [T^) with that of its “multi-SSP” set of simula¬ 
tions. The overall total number of stars in both sets of SSP sim¬ 
ulations is normalised to the number of cluster stars on the CMD 
brighter than the 50% completeness limit. Finally, all individual 
simulations are repeated 20 times and their cross-SGB profiles are 
averaged together. 

3.2 Results 

3.2.1 NGC 1651 

The resulting cross-SGB profiles for NGC 1651 are shown in Fig¬ 
ure HJi. The profiles were derived as density functions, using the 


non-parametric Epanechnikov kernel iSilvermanll986h . This mini¬ 
mizes biases that can arise if fixed bin widths are used. 

The cross-SGB profile for the cluster itself (black solid line) 
shows some significant differences with that of the single-SSP sim¬ 
ulation (blue line): the cross-SGB profile of the cluster is signifi¬ 
cantly wider than that of the single-SSP simulation, and there is an 
offset between the two profiles in that the cluster profile is centered 
on a F475W magnitude that is ~ 0.08 mag brighter than that of the 
single-SSP simulation. The same offset is present between the pro¬ 
files of the cluster and that of its “multi-SSP” simulation (see red 
solid line), while the width and shape of the latter are consistent 
with that of the cluster’s profile (see dashed red line, which is off¬ 
set by 0.08 mag from the solid red line in the positive X direction). 
The similarity of the shapes of the cross-SGB profiles of the cluster 
and that of its multi-SSP simulation, along with the profile of the 
single-SSP simulation being significantly narrower, strongly sug¬ 
gests that the SGB morphology of NGC 1651 is better described by 
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Figure 3. Panel (a): The black solid line represents the cross-SGB profile 
for NGC 1651, while the light grey region indicates its 68% confidence 
interval. The red solid line represents the cross-SGB profile of the “multi- 
SSP” simulation described in Sect. EH 1. The red dashed line is the same 
after shifting it by 0.08 mag to the right. For comparison, the solid blue line 
represents the cross-SGB profile of a single-SSP simulation of NGC 1651, 
divided by a factor of ~ 3.5 in order to scale its maximum density to that 
of NGC 1651. Panel (by. Comparison between “multi-SSP” simulations of 
NGC 1651 created with binary fractions 0.0 (long dashed line), 0.2 (solid 
line), and 0.4 (short dashed line). 


a distribution of ages similar to that implied by its cross-MSTO pro¬ 
file, rather than by a single SSP. This conclusion is opposite to that 
of L-l-14, who stated that the cross-SGB profile of NGC 1651 “can 
be explained only by a single-age stellar population, even though 
the cluster has a clearly extended main-sequence turn-off region” 
(quoting their abstract). Since both conclusions were drawn from 
the same original dataset, we discuss this disagreement in some de¬ 
tail below. 

The analysis methods of L-l-14 differ in two main ways from 
those of ours. One difference is that L-l-14 did not use SSP sim¬ 
ulations in their analysis. The ages in their cross-SGB profiles 
(their Fig. 4) were determined directly from the location of the 
isochrones in the CMD, thus only involving single stars. Con¬ 
versely, we use SSP simulations including measurement errors and 
unresolved binary stars in this context. We note that the binary 
fractions of inter mediate-age star clu s ters like NGC 1651 are sub¬ 
stantial (» 20%. I^kev et all [200^: iMilone et al.ll2009l: G-flla; 
G-l-14), and our simulations show that binaries extend the luminos¬ 


ity distribution of the SGB towards brighter magnitudes at the ages 
of such clusters (see, e.g.. Fig. [Tjt)- The effect of the binary frac¬ 
tion on the cross-SGB profile shape is illustrated in Fig.[^ for bi¬ 
nary fractions of 0%, 20%, and 40%. Secondly and perhaps more 
importantly, L-l-14 did not take the stars in the SGB region of the 
CMD located above their log [age (yr)] = 9.24 isochrone into ac¬ 
count (see their Figs. 2 and 4). Given our results shown in Fig. [3^, 
we suspect that if they would have done so, they would have found 
that the cross-SGB profile of NGC 1651 continues up to brighter 
magnitudes (equivalent to younger ages) with a profile shape that 
is clearly wider than that of a single SSP. 

Finally, the offset in magnitude between the cross-SGB profile 
of NGC 1651 and that of its SSP simulations is also likely rel¬ 
evant to the disagreement between the conclusions of L-l-14 and 
ours, since L-l-14 based their conclusions in part on the observation 
that the age distribution measured from the cross-SGB profile of 
NGC 1651 reached an age lower than that measured from its cross- 
MSTO profile (cf. their Fig. 4). This will be discussed below. 

3.2.2 IC2146 

We recall that this cluster was included in our sample because its 
mass and age are similar to those of NGC 1651; while NGC 1651 
features an eMSTO, IC 2146 does nojj As such it provides a rele¬ 
vant comparison in the context of the nature of the SGB morpholo¬ 
gies of such clusters. 

Figure |4jj shows a comparison between the cross-SGB pro¬ 
files of IC 2146 (black line) and that of its single-SSP simulation 
(blue solid line). Note that in this case, the width of the profile of 
the single-SSP simulation is very similar to that of the cluster it¬ 
self. This difference with respect to NGC 1651 is consistent with 
IC 2146 not hosting an eMSTO whereas NGC 1651 does, although 
this may in principle also be due in part to the larger photomet¬ 
ric errors in the dataset of IC 2146 relative to those of the other 
clusters presented in this paper. However, an important result of 
the simulations of IC 2146 is that the magnitude offset between the 
cross-SGB profiles of NGC 1651 and its SSP simulation (cf. above) 
is also present for 1C 2146. This suggests that the cause of this off¬ 
set is not related to the question whether or not an intermediate-age 
cluster hosts an eMSTO (or a distribution of stellar ages). Instead, 
we suggest that the offset is of a systematic nature, related to one 
or more parameters relevant to the MSTO-SGB era of stellar evo¬ 
lution theory for stars of moderate metallicity and masses of « 1.5 
M0. This is explored further in Sect. [4] 

3.2.3 NGC 1783, NGC 1806, and NGC 1846 

We performed a similar analysis of the SGB regions in the massive 
eMSTO cl usters NGC 1783 NGC 1806, and NGC 1846, using the 
data from iGoudfrooii et al.l l l2009l) and G-l-lla. Since the ages of 
these three clusters are very similar to one another, we adopt the 
M-l-08 isochrone for an age of 1.65 Gyr to set the magnitude ze- 
ropoint for the cross-SGB profiles of these clusters. The results 
are shown in Figs. 0-13 Similar to the case of NGC 1651, we 
find that the shapes of the cross-SGB profiles of all three clusters 
are consistent with those of their respective multi-SSP simulations, 
while they are significantly wider than those of their single-SSP 
simulations. And again, the clusters’ cross-SGB profiles are offset 

^ G-l-14 suggested that this difference is due to their different escape veloc¬ 
ities. 
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Figure 4. Panel (a): CMD of IC 2146 within its core radius, taken from 
ICorrenti et alJ i2Q14h . along with a M+08 isochrone of age 1.90 Gyr. The 
SGB selection box is outlined by black dashed lines. Red dots indicate stars 
in the outskirts of the HST image with a total area equal to that within 
the core radius of IC 2146. Panel (by. The black solid line represents the 
cross-SGB profile for IC 2146, while the light grey region indicates its 68% 
confidence interval. For comparison, the solid blue line represents the cross- 
SGB profile of a single-SSP simulation of IC 2146. The blue dashed line is 
the same after shifting it by 0.13 mag to the right. 


Figure 5. Panel (a)\ CMD of NGC 1783 within its core radius, taken from 
G+1 la, along with M-l-08 isochrones of age 1.50, 1.60, 1.70, 1.80, and 1.90 
Gyr (top to bottom). Red dots indicate stars in the outskirts of the HST 
image with a total area equal to that within the core radius of NGC 1783. 
The SGB selection box is outlined by black dashed lines. Panel (b): Similar 
to Figure[3] but now for NGC 1783. The SGB profile of the “Single SSP 
Simulation” was divided by a factor of ~ 4.5 in order to scale its maximum 
density to that of NGC 1783. 


from those of their multi-SSP simulations by a shift of 0.08 mag. 
Similar to the interpretation of this offset by L+14 for the case of 
NGC 1651, BN15 interpreted this offset for the cases ofNGC 1806 
and NGC 1846 in the sense that their SGB morphologies are incon¬ 
sistent with the presence of a significant age spreaqj. 


Note however that the shapes of the SGB profiles ofNGC 1806 and NGC 
1846 shown by BN15 (see their Figs. 3 and 8) seem to be more consistent 
with their SSP simulations that involve a spread of ages than with those of 
a single-age SSP. 


3.2.4 Summary 

Summarizing the results in this Section, it seems that the shapes of 
cross-SGB profiles of eMSTO clusters are consistent with those of 
their respective multi-SSP simulations, while they are inconsistent 
with those of their single-SSP simulations. Note that in this sense, 
the cross-SGB profiles of such clusters are entirely consistent with 
their cross-MSTO profiles when the latter are interpreted as a distri¬ 
bution of ages (i.e., the pseudo-age distributions derived by G+1 la 
and G14). 

The main inconsistency between the cross-SGB profiles of 
eMSTO clusters and those of their respective multi-SSP simula¬ 
tions shown above is that of an magnitude offset of 0.1 mag be¬ 
tween them (for the filters of the HST instruments used here, i.e., 
F475W of WFC3 and F435W and F555W of ACSAVFC) in that 
the clusters’ SGB stars are brighter than the simulations that use 
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Figure 6. Similai- to Figure|5] but now for NGC 1806. The SGB profile of 
the “Single SSP Simulation” in panel (b) was divided by a factor of ~ 4.6 
in order to scale its maximum density to that of NGC 1806. 



Figure 7. Similar to Figure[5] but now for NGC 1846. The SGB profile of 
the “Single SSP Simulation” in panel (b) was divided by a factor of ~ 4.4 
in order to scale its maximum density to that of NGC 1846. 


their best-fitting M+08 isochrones. This issue is addressed in the 
next Section. 


4 DISCUSSION 

4.1 The Impact of Convective Overshoot 


Stars in the MSTO and SGB region of the CMD for ages of ~ 1.5 - 
2.5 Gyr have masses in the range 1.3 ^ (M/Mq) ^ 1.6. One pa¬ 
rameter in stellar evolution theory that is known to have a signifi¬ 
cant impact on the photometric properties of stars in this mass range 
is that of the level of “overshoot” from the convective core, i.e., 
the mean free path of c onvective bubbles across the border of the 
convective region (e.g., Meade^ 19751: Bressan. Chiosi. & Bertellil 


Il98ll : iBertelli. Bressan. & Chiosil 1 19851 : IPemargue et al.l l2004h . 
The level of convective overshoot Ac depends on stellar mass in the 
mass range between fully radiative cores {M SS 1.0 Mq) and that 
of fully convective cor es (M ^1.5 Mp)), and it is also believed t o 
depend on metallicity dPemaraue et al] 20041 : iBressan et al ]|2 oT^ . 
The detailed dependencies of Ac on stellar mass and metallicity 


are not well-established, due in part to a relative lack of high- 
quality data of massive intermediate-age star clusters for calibra¬ 
tion purposes. A consequence of this uncertainty is that different 
SSP models implement the dependencies o f Ac on stellar mass and 
metallicity in different ways (compare , e.g., Pietrinferni_et alj2004 


IPotter et^l2008l: iMarigo et al.ll2008l: IBressan et al.ll2012 ). 

In the following we investigate whether an adjustment of the 
level of overshoot can reasonably account for the observed mag¬ 
nitude offset between the SGB profiles of clusters and those of 
their SSP simulatio ns. To do so, we run a set of isochrones based 
on PARSEC tracks dBressan et alj|20l3l using a range of Ac. The 
model set covers 0.3 0 < Ac < 0.70, where values for Ac are in the 
IBressan et al.l l ll98lh formalisnj^, for a range of ages and metallic- 
ities relevant to the clusters presented in this paper. Note that there 
are additional ways to adjust the level of overshooting for a given 
isochrone, such as shifting the stellar mass range in which over- 


^ Ac refers to the maximum level of overshoot in the relevant stellar mass 
range. The increase of overshoot level with stellar mass is assumed to be 
linear. 
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shooting tak es place or by adop t ing a diffusiv e approach for this 
process fe.g.. lFrevtag et al.lll99d : lHerwidl200(]li . 

Fig. [8^ shows the CMD of NGC 1651 along with three PAR¬ 
SEC isochrones with levels of overshoot ranging from Ac = 0.35 to 
0.50. For each case, the age of the isochrone was adjusted to coin¬ 
cide with the best-fitting isochrone from the M-l-08 models (shown 
in Fig. 12 and with a black solid line in Fig. HI near the red end of 
the turnoff region (at F475W ~ 21.15 and F475W-F814W ~ 0.82 
to be precise). The values for (m-M)o and Ay were chosen to re¬ 
produce the positions of the RGB, the lower MS, and the red clump 
(RC), in a way equivalent to that described by G-rlla. Note that 
the shape and location of the top of the MSTO of the PARSEC 
isochrone with Ac = 0.40 are virtually identical to those of the best¬ 
fitting M-r08 isochrone, while its SGB is clearly brighter (by up to 
0.12 mag at the red end of the SGB). Apart from this difference at 
the SGB and a small colour difference at the bottom of the MSTO, 
the PARSEC isochrone with Ac = 0.40 provides a fit to the CMD 
features of NGC 1651 that is very similar to that of the best-fitting 
M-l-08 isochrone. Fig.[2> illustrates the impact of the use of this 
PARSEC isochrone with Ac = 0.40 for creating cross-SGB profiles 
of NGC 1651 and its SSP simulations. This plot represents a copy 
of Fig.[2except that the zeropoint of the SGB magnitude offsets is 
now set by the PARSEC isochrone rather than the best-fitting M-l-08 
isochrone. Note that the cross-SGB profile of NGC 1651 is now 
fully consistent with that of its multi-SSP simulation (i.e., without 
any residual shift), while its single-SSP simulation is still clearly 
narrower. 

Given this result, we suggest that the observed magnitude off¬ 
sets between the SGB profiles of star clusters and those of their SSP 
model isochrones or simulations can be understood entirely by de¬ 
tails of the way the dependencies of Ac on stellar mass and metal- 
licity are implemented in the SSP model being used. As such, we 
argue against the use of the presence of such (small) magnitude off¬ 
sets between cross-SGB profiles to confirm or deny the presence of 
a given age distribution in intermediate-age star clusters such as the 
ones discussed in L-l-14, BN15, and the current paper. Instead, we 
suggest that comparisons between SGB properties of such clusters 
and their SSP model predictions be focused on the shapes of their 
cross-SGB profiles, for clusters with statistically sufficient numbers 
of stars on the SGB, and that unresolved binaries are properly taken 
into account in such analyses. 


4.2 The width of the Red Clump as a diagnostic for age 
spreads 

In addition to their study of the shapes of the SGB, BN 15 also pre¬ 
sented a comparison between the RCs in NGC 1806 and NGC 1846 
and those of M-l-08 isochrones of different ages (their Pigs. 4 and 9), 
concluding (1) that the observed RC’s showed a smaller spread in 
colour than expected from clusters having a significant age spread, 
and (2) that the isochrones that fit the youngest part of the eMSTO 
were also those that showed the best fit to the mean RC position. 

To put BN15’s analysis in context, we call attention to a few 
points that are relevant to the interpretation of the RC’s position 
and width. Pirst, the (mean) position of the RC depends on a num¬ 
ber of factors in addition to age, distance, and reddening. One such 
factor is the mixing-length parameter qmlt, which defines the ef¬ 
ficiency of energy transport in the external layers of the star and 
hence its effective temperature Teff for a given luminosity. In most 
isochrone sets, including the M-l-08 ones, aMLT is calibrated in the 
solar model, and then applied to all stars including the r ed gia nts. 
Recent 3-D model atmospheres from iTrampedach et al.l l l2014h in- 
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Figure 8. Panel (a)\ CMD of NGC 1651 showing the upper MS, MSTO, 
SGB, lowe r RGB, and Red Cl ump features. The black line represents the 
best-fitting lMarigo et al.l il^OOSl) isochrone (cf. Fig.[^). The green, red, and 
blue lines represent PARSEC isochrones with [M/H] = —0.40 for which 
the ages and overshoot values are indicated in the legend. Panel (b): Same 
as Fig. [5] but now using the SGB of the PARSEC isochrone with age = 
2.2 Gyr and Ac = 0.40 as the zeropoint of the cross-SGB profiles. Note 
the similai'ity of the cross-SGB profiles of NGC 1651 (black line with 68% 
confidence interval in light grey) and that of its “multi-SSP simulation” (red 
line), while the cross-SGB profile of the “single-SSP simulation” (blue line) 
is much narrower. 


dicate however that while the approximation of a nearly constant 
ctMLT might be good for stars along the upper RGB (and probably 
also for the RC), it also indicates significant variations between the 
ctMLT values of red giants and those of dwarfs such as those popu¬ 
lating the MS of Magellanic Cloud clusters. The variations are such 
(a few tenths in ctMLi) that they could change the relative colours 
of MS and RC stars by a few hundredths of a mag. Another uncer¬ 
tainty in isochrone colours is due to the adopted Teff-colour rela¬ 
tions, which_arejasuany(de^ from 1-D static model atmospheres 
(e.g. JCastelli & Kuruczll2Q0^ . Therefore, one should be cautious 
in interpreting issues with regard to ages of star clusters based on 
differences between the MS and RC positions in the CMD, relative 
to isochrones, at levels of 0.01-0.03 mag. 
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As to the width of the observed RC, we agree with BN15 that 
it is compact in both NGC 1806 and NGC 1846, with its main com¬ 
ponent spanning just « 0.5 mag in F435W-F814W colour. FIow- 
ever, both NGC 1806 and NGC 1846 also exhibit a number of sec¬ 
ondary RC stars (hereafter SRC: lGirardi et ^l2009l : G-l-14). These 
are He-burning stars that were massive enough to avoid electron 
degeneracy settling in their cores after leaving the main sequence. 
In both clusters, the SRC stars define an almost vertical feature to 
the bottom left of the main RC, centered at F435W-F814W ~ 1.80 
mag, and with F814W magnitudes between 18.4 and 18.7 (see Fig. 

they are seen also in BN15’s plots and are somewhat better de¬ 
lineated in NGC 1806 than in NGC 1846. These SRC sequences 
are, by themselves, a clear indication that the RC in these clus¬ 
ters is not as homogeneous and compact as claimed by BN15. In 
the age-spread interpretation of eMSTOs, the presence of a SRC 
is simply reflecting the spread in ages (i.e., initial stellar masses) 
of stars leaving the MSTO, and indeed the fraction of RC stars in 
the SRC among eMSTO clusters correlates strongly with the frac¬ 
tion of MSTO stars at the youngest (brightest) end of the MSTO 
(G-l-14). In the rotating scenario, instead, SRC stars are the putative 
progeny of fast MS rotators, since rapidly rotating stars have larger 
core masses at the end of the MS era than do nonrot ating stars (e.g., 
iMaeder & Mevnedbood : l^genberger et alJbOlOl) . However, this 
does not seem to explain the aforementioned observed correlation 
with the brighter part of the eMSTO (see G-l-14 for more details). 

That said, it is interesting to note that such SRCs in NGC 
1806 and NGC 1846 apparently contrast with comparatively nar¬ 
row “main” RCs. As noticed by BN15, M-l-08 isochrones of ages 
between 1.3 and 1.8 Gyr span a significant range of RC colours 
(A(F435W-F814W) ~ 0.1 mag), which seems wider than the RC 
colour spreads in their Figs. 5 and 10. There are however two ad¬ 
ditional points to be taken in consideration in this context: firstly, 
BN15 inferred the age distribution from the stars in the RC by as¬ 
signing ages based on the “closest isochrone” rather than using a 
smooth function. This method introduces an intrinsic binning in 
the age distribution, with a bin size depending on the age spac¬ 
ing of the isochrones. Secondly, the M-l-08 isochrones were created 
by the interpolatio n of evolutionary tracks (originally taken from 
iGirardi et alj|200d) which, for the mass interval into consideration 
(A4 ~ 1.5- 1.8 Mq), were computed for initial masses spaced by 
AM = 0.05 Mq. The issue with such isochrones in the context of 
the interpretation of the RC morphologies of these star clusters is 
that they cannot reveal any CMD structure caused by variations 
in stellar evolutionary behaviour that occur over a mass interval 
smaller than 0.05 Mq. If any such variations occur, the isochrones’ 
behaviour as a function of age will reflect the smooth interpola¬ 
tion between the confining tracks, rather than the real expected 
behaviour. This is what we believe is happening with the M-l-08 
isochrones. 

As recently noticed by IGirardi et al.l l l2013bl . see their ap¬ 
pendix), the onset of electron degeneracy in stellar cores at A4 
1.7 Mq (with the exact limiting value depending on overshooting 
efficiency) produces a marked discontinuity in the evolutionary be¬ 
haviour of stars, occurring over a mass scale smaller than 0.01 Mq. 
This significantly affects the expected shape of the RC in more de¬ 
tailed isochrone models, where stars are in principle expected to be 
either on the main RC or in the SRC, but rarely in between these 
features (modulo unresolved binary stars). As a consequence, the 
RC is expected to be more compact than that displayed in M-l-08 
isochrones. The expectation from PARSEC isochrone models com¬ 
puted from tracks with a fine grid of stellar masses - namely those 
derived from a Z = 0.008 grid with the same level of detail and 


mass resolution as in IGirardi et al.1 l l2013hl . i.e., AM/Mq = 0.01) 
- is shown in Fig. |9ji, and compared with M-l-08 isochrones (Fig. 
Igjt) using the same values for (m-M)o and A^. In both cases, we su¬ 
perimpose the stars observed within the core radius of NGC 1806. 
It can be seen that, with the new PARSEC models, the extent of the 
observed “main” RC appears to be compatible with an age range 
of Alog(age) ~ 0.10 (or A(age) ~ 500 Myr at an average age of 
1.7 Gyr), whereas it indicates a smaller age range in the case of the 
M-l-08 isochrones (Alog(age) ~ 0.07 or A(age) ~ 200 Myr). 

The ability of the new PARSEC isochrone models that were 
computed from a grid of initial stellar masses spaced by AM = 
0.01 Mq to fit the full CMD of NGC 1806 is illustrated in Fig. 
13;. We used (m-M)o = 18.53 mag and Av = 0.035 mag for these 
isochrones. Note the good fit to the MSTO, SGB, RGB, and RC 
regions. 

Finally, we determine the age distribution from the RC of 
NGC 1806 by means of the distribution of stars in a rectangle in 
the RC region of the F814W vs. F435W-F814W CMD. The short 
and long axes of the rectangle are approximately parallel and per¬ 
pendicular to the isochrones, respectively (see blue rectangle in Fig. 
[3>). The magnitudes and colours of stars in the rectangle are then 
transformed into the coordinate frame defined by its two axes. The 
star positions in the coordinate perpendicular to the isochrones are 
then translated to age by repeating the same procedure for the new 
PARSEC isochrone tables described above, for an age range that 
fully covers the observed extent of the RC region of NGC 1806 (us¬ 
ing the same values of Z, (m-M)o and Ay), and performing a poly¬ 
nomial least-squares fit between age and the coordinate perpendic¬ 
ular to the isochronefl The resulting age distribution of NGC 1806 
is shown in Fig.[3i (black line), using a non-parametric Epanech- 
nikov density kernel as before. For comparison, Fig.[3i also shows 
the age distribution that would be expected for a single SSP, derived 
from a set of Monte Carlo simulations (as described in Sect. I3.lt 
using the new PARSEC isochrone for an age of 1.74 Gyr. Note that 
while the age distribution is shifted to slightly older ages relative 
to those indicated by the cross-MSTO profiles when the latter are 
compared with M-l-08 isochrones (cf. Fig.[T};), the cross-RC profile 
shape of NGC 1806 is similar to that of its cross-MSTO one, and it 
is significantly wider than that of a single-SSP simulation. 

Although these new RC models will be described in more de¬ 
tail in a forthcoming paper, these results already demonstrate that 
the expected detailed RC morphology of intermediate-age clusters 
that host a range of stellar ages is model -and mass-resolution- de¬ 
pendent. We therefore argue that the conclusions by BN15 on the 
nature of the RCs of NGC 1806 and NGC 1846 were largely due 
to their choice of isochrones and analysis methods. In contrast, our 
results with newly updated isochrone models strongly suggest that 
the RC morphologies of these clusters are actually described very 
well by a distribution of stellar ages that is consistent with that 
indicated by their MSTOs. 

4.3 Implications regarding the nature of eMSTOs 

Our result that the widths and shapes of cross-SGB profiles of eM¬ 
STO clusters are consistent with those of their respective multi-SSP 

® We excluded the SRC area of the CMD from this rectangle to avoid the 
highly non-linear (and multi-valued) functional relation between age and 
the coordinate along the long axis of such a rectangle in the SRC region. 
This causes an underestimate of the number of stars at the youngest ages 
shown in Fig. HJl. 
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Figure 9. Panel (a): The CMD of NGC 1806 (cf. Fig.|^), zoomed in on the RC region. The lines renresent lMarigo et al.l i2008h isochrones for log (age/yr) = 
9.10 to 9.24 (left to right) with a spacing of 0.02. Panel (b): Same as panel (a), but the lines now represent new PARSEC isochrones with Z = 0.008 which were 
derived from stellar tracks spaced by AA1 /Mq = 0.01 (see Sect. 14.2t . Isochrones are for log (age/yr) = 9.18 to 9.32 (left to right) with a spacing of 0.02. The 
thick blue rectangle in panels (a) and (b) is used to derive the age distribution of RC stars as described in Sect. EJ Panel (c): CMD of NGC 1806 showing the 
region from below the MSTO to beyond the RC along with new PARSEC isochrones for Z = 0.008 and ages of 1.58, 1.74, 1.91, and 2.09 Gyr (left to right). 
Panel (d): The black solid line shows the cross-RC distribution of NGC 1806 expressed in age units (see Sect. 112), while the light grey region indicates its 
68% confidence interval. For comparison, the solid blue line represents the cross-RC distribution of a single-SSP simulation of NGC 1806, divided by a factor 
~ 1.9 to yield the same maximum star density as that of NGC 1806 itself. 


simulations, and significantly wider than those of their single-SSP 
simulations, has a number of interesting implications related to the 
nature of eMSTOs. 

It is obvious that this result is consistent with the prediction 
of the “extended star formation” scenario; we therefore focus on 
a comparison with predictions of the “stellar rotation” scenario. 
It is well known that the centrifugal force in rotating stars with 
masses in the approximate range of 1.2-1.7 Mq decreases their 
effec tive gravity which in tur n lowers their e ffective temperature 
(e.g., Mevnet & Maedeilll99'^ . This prompted ISastian & de Minlj 
1 20091) to suggest that a range of stellar rotation velo cities may 
cause the eMSTO feature. However. iGirardi et al.l ( l201lh computed 
an isochrone for a typical age of eMSTO clusters and a typical stel¬ 
lar rotation velocity and found that the prolonged lifetime of rotat¬ 
ing stars approximately cancels the effects of the lower gravity in 
terms of the MSTO colour, which was found to be virtually equal 
to that of the non-rotating isochrone. The main impact of stellar 


rotation was found to be a lengthening of the MSTO hook in con¬ 
juncti on with a brightening of the SGB (see Fig. 3 in iGirardi et al.l 
l201lh . similar to the effects of an increased level of convective 
overshoot shown in Sect. ML. _. 

More recentlv. fyang et al] ( l2013h computed isochrones that in¬ 
cluded the effects of stellar rotation, using the Yale Rotating Evolu¬ 
tion Code dPinsonneault et alii 19891) . They found that rotation can 
actually cause eMSTOs if the efficiency of rotational mixing is suf¬ 
ficiently smalfl. However, the luminosities of the onset of the SGB 
of the rotating isochrones was found to be consistent with those of 
the corresponding non-rotating isochrones. This wa s also part of 
the arguments of L-l-14, who used the recent tracks of lGeorgy et al.l 


^ Note however that G-l-14 pointed out that the width of the MSTO of the 
most massive eM STO clusters is significantly larger than the predictions of 
lYangetalJ <20131) . even for their smallest value of rotational mixing effi¬ 
ciency. 
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(l2014h and showed that the track for a rapidly rotating star is at a 
cooler temperature than that of the non-rotating track at the MSTO, 
while the rotating track merges with the non-rotating track right 
at the onset of the SGB. However, we remind the reader that due 
to the longer lifetime of rotating stars, one should not use stellar 
tracks to interpret the morphologies of the MSTO or SGB in CMDs. 
One should compute isochrones instead. This will be addressed in 
a forthcoming paper. 

In summary, a spread of rotational velocities in a single-age 
stellar population would cause the MSTO and SGB features to 
show the following properties according to the predictions of the 
recent modeling efforts: 

• According to the findings of iGirardi et alj (1201 ih . the MSTO 

would not show any significant spread in colour, only a spread 
of the length of the MSTO hook. This is inconsistent with the 
MSTO morphologies of eMSTO clusters. Conversely, the SGB 
would show a significant spread of luminosities (similar to what 
we find in eMSTO clusters). _ 

• According to I Yang et al.l j2013h and L-l-14, the MSTO would 
show some significant spread in colour (similar to the morphology 
of eMSTOs) if the efficiency of rotational mixing is small, but the 
luminosity of the SGB would not be affected in a significant way. 
The latter prediction is inconsistent with our cross-SGB profiles of 
eMSTO clusters. 

Given the above, our result that the cross-SGB profiles of eM¬ 
STO clusters are significantly wider than that of their single-SSP 
simulations (and consistent with the widths of MSTOs of eMSTO 
clusters) seems to indicate that their MSTO and SGB morpholo¬ 
gies are better described by a distribution of stellar ages than by 
the combination of a single age and a range of stellar rotation ve¬ 
locities. We recognize that the extent to which rapid stellar rota¬ 
tion affects the temperature and luminosity of SGB stars is still in 
some l evel of flux. For example, the grid of models bv ICeorgy et al.l 
( I2OI4I) has so far only been completed for stars with Ad > 1.7 Mq, 
while SGB stars in the eMSTO clusters discussed in this paper are 
less massive (by 0.1 -0.3 Mq As such, our conclusions refer to 
the situation at the time this paper is written. 


5 SUMMARY AND CONCLUSIONS 

In the context of the question of the nature of eMSTOs in massive 
intermediate-ag e star clusters in the LM C, we investigated the 
recent claims of lUi et aP ( 1201 4 L-(-14) and iBastian & NiederhofeJ 
( I2OI5L BN15) regarding the nature of such clusters’ SGB mor¬ 
phologies, which are expected to be affected more by age spreads 
than by spreads in stellar rotation velocity. Their analysis led them 
to argue that the SGB morphologies of three clusters featuring 
eMSTOs (NGC 1651, NGC 1806, and NGC 1846) are inconsistent 
with extended star formation histories within those clusters. We 
performed an independent study of the SGB morphologies of five 
intermediate-age star clusters in the LMC, including the three 
clusters studied by L-(-14 and BN15 (NGC 1651, NGC 1806, and 
NGC 1846). Comparisons between the SGB morphologies of 
cluster stars and those of SSP model predictions are done using 
Monte-Carlo simulations of single-age SSPs as well as “composite 
SSPs” whose age distributions are taken from the “pseudo-age 
distributions” of the clusters which were derived from their MSTO 

* We note that the study of I Yang et al.l fe013l) showed that the impact of 
stellar rotation to isochrone shapes changes significantly in this mass range. 


morphologies by G-i-lla and G-(-14. The SSP simulations include 
a self-consistent treatment of unresolved binaries. Other method¬ 
ological differences with the L-(-14 and BN15 studies are pointed 
out where relevant. Our main conclusions are the following. 

(i) In contrast with L-(-14 and BN15, we find that the shapes of 
the cross-SGB profiles of all clusters in our sample are consistent 
with their cross-MSTO profiles (to within la) when the latter are 
interpreted as distributions in age. Conversely, cross-SGB profiles 
of simulated single-age stellar populations (convolved with pho¬ 
tometric uncertainties) are found to be significantly narrower than 
those of the clusters. 

(ii) We argue that the magnitude offset between the cross-SGB 
distributions of the cluste rs in our sample an d that of their best¬ 
fitting isochrones from the lMarigo et al.l ( l2008h family of SSP mod¬ 
els can be understood once the variation of the treatment of convec¬ 
tive overshoot among SSP isochrone models is taken into account. 

(iii) We investigate the claims made by BN 15 who stated that 
the red clumps of NGC 1806 and NGC 1846 showed a smaller 
spread in colour than expected from clusters having a significant 
age spread, and that the M-(-08 isochrones that fit the youngest part 
of the eMSTO were those that showed the best fit to the mean RC 
position. We find that these results are dependent on the models and 
methods used. Using isochrones newly created from a grid of PAR¬ 
SEC tracks that features a 5-fold higher resolution in stellar mass 
than that used by the M-f08 (and most other sets of) isochrones, 
we find that ages indicated by the MSTO of NGC 1806 are actu¬ 
ally consistent with those indicated by the RC. Furthermore, we 
find that the RC morphology of NGC 1806 is consistent with the 
age distribution indicated by the MSTO, and significantly more ex¬ 
tended than that expected from a single-age SSP. 

(iv) We compare the observed MSTO and SGB morphologies 
of eMSTO clusters with those of (current) predictions of models 
that assume (a) a distribution of stellar ages versus (b) a range of 
stellar rotation velocities at a single age. This comparison indicates 
that a distribution of stellar ages is currently the only option that 
can explain both the observed MSTO and SGB morphologies. 

Our overall conclusion is that in spite of recent arguments by L-l-14 
and BN 15, the SGB and RC morphologies of star clusters featuring 
eMSTOs are consistent with the scenario in which eMSTOs are 
(mainly) due to a distribution of stellar ages. 
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